Gauge fields at the surface of topological insulators 



M. I. Katsnelson 
Radhoud University Nijmegen, Institute for Molecules and Materials, 
Heyendaalseweg 135, NL-6525 AJ Nijmegen, The Netherlands 

F. Guinea and M. A. H. Vozmediano 
Instituto de Ciencia de Materiales de Madrid (CSIC), Cantoblanco, Madrid 28049, Spain 

(Dated: January 26, 2013) 

We study the emergence of gauge couplings in the surface states of topological insulators. We 
show that gauge fields arise when a three dimensional strong topological insulator is coated with 
an easy-plane ferromagnet with magnetization parallel to the surface. We analyze the modification 
induced by the gauge fields on the surface spectrum for some specific magnetic configurations. 



PACS numbers: 

INTRODUCTION 

The Dirac fermion description of novel condensed mat- 
ter systems have given rise to a great deal of new phe- 
nomena and ideas in the field that is acting as a bridge 
between different areas of physics including quantum field 
theory, cosmology, elasticity, and statistical mechanics. 

One of the most successful examples can be found in 
the physics of graphene [T] and, in particular, in the at- 
tractive relation between the structural - lattice - prop- 
erties and the electronics of this material. It is by now 
well known that elastic deformations and curvature of 
the graphene lattice give rise to effective gauge fields that 
couple to the fermionic degrees of freedom with the min- 
imal coupling typical of gauge theories [5] . 

Another important example of Dirac fermions is that of 
topological insulators. These materials are characterized 
by a bulk insulating behavior with conducting bound- 
ary states protected by some topological property, the 
integer quantum Hall effect |S1 Sj being the prototype. 
In the recent developments (for a review, see Refs. \^ 
[7]) numerous systems have been found where the edge 
states exist in the absence of a magnetic field and are 
protected by time reversal symmetry. The simplest ex- 
ample is the spin Hall effect [51 [S] where the spin orbit 
coupling plays a major role; most of the newly discovered 
topological insulators are based on materials with strong 
spin orbit coupling [T0lfl2] . The topological protection of 
edge states ensures their stability against non-magnetic 
- more generally, time-reversal invariant - perturbations, 
opening up the possibility for future spintronic devices. 

Novel exotic states appear at the surface of a three 
dimensional (3D) topological insulator when an energy 
gap is induced by a perturbation breaking the time re- 
versal symmetry. They are obtained typically by apply- 
ing a magnetic field perpendicular to the surface what 
gives rise to the half integer quantum Hall effect [TSj or 
by proximity effects with a magnetic material what in- 
duces an anomalous quantum Hall effect [Mj . A very re- 
cent experiment reports the observation of massive edge 



Dirac fermions in a three dimensional topological insula- 
tor doped with magnetic impurities |15j . In these cases 
the exotic physics is due to a combination of the regular 
or space-dependent mass induced by the perpendicular 
magnetic field to the Dirac fermions, and by the appear- 
ance of zero modes associated to the Dirac operator in 
magnetic fields [I6l418j . 

A different situation arises when the surface of the 
topological insulator is coated with an insulating ferro- 
magnet with an easy plane such that the magnetization is 
parallel to the surface. Although time reversal invariance 
is also broken and edge states are no more protected, a 
mass is not directly induced in the system. In this work 
we will see that in such situation gauge fields appear cou- 
pling to the electronic degrees of freedom similar to these 
produced by elastic deformation on graphene. We will 
explore the modification of the spectrum for the cases of 
a vortex and a skyrmion configuration. 



GAUGE FIELD INDUCED BY AN IN-PLANE 
MAGNETIZATION 

When putting a thin ferromagnetic layer with magne- 
tization M on the surface of a strong topological insula- 
tor an interaction term of the form H' = —JM.a arises 
where J is the coupling constant. There are two contri- 
butions to this magnetic term: the Zeeman energy, and 
s—d exchange interaction energy. The latter contribution 
can be orders of magnitude higher than the former one 
if the overlap between the electron wave function on the 
surface of the topological insulator and that of the ferro- 
magnet is big enough. One can easily expect that JM ^ 
0.01-0.1 eV. This term breaks the time reversal symme- 
try. In the case of an easy-axis ferromagnet with the 
magnetization perpendicular to the layer the magnetic 
interaction has a structure of a mass term proportional 
to az in the Dirac Hamiltonian resulting in the gap open- 
ing mentioned above. In the case of having an easy-plane 
ferromagnet, the vector M lies in the surface plane and 
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FIG. 1: (Color online) Sketch of the ferromagnetic vortex 
configuration. 



the s — d exchange interaction couples to the fermions as 
a vector potential. 

In the case of graphene where the "spin" degree of 
freedom is actually a pseudospin related to the sublat- 
tice degree of freedom, gauge couplings are induced by 
strains or other mechanical deformations of the lattice 
[2]. In topological insulators where spin is the real spin 
one needs in-plane magnetization to create such terms. 

Easy plane fcrromagnetism is very natural for thin 
films of single layers as the shape anisotropy (magnetic 
dipole-dipole interaction) always favor this situation |19j . 
This will be the case for any soft magnetic material i. e. 
materials fulfilling the condition AttM"^ > K, K being the 
constant of magnetocrystalline anisotropy. If the coating 
is a circular symmetric disc, the magnetization vector will 
be freely rotating in the plane. Magnetic textures at the 
surface of a three dimensional topological insulator have 
been recently studied in Refs. 20, 21, In the latter pub- 
lication, the exchange field energy associated to parallel 
ferromagnetic domains of 50 nm width was estimated to 
be of the order of 6 meV. 



Vortex configuration 

In easy-plane ferromagnets vortices can exist which are 
spin configurations with a net 2n twist about a particular 
point or core similar to the Onsager-Feynman vortices in 
superfluid helium or Abrikosov vortices in superconduc- 
tors. Examples are BaCoAs04 and K2CUF4 22 . 

In this case the magnetization is constant in magnitude 
and rotates in direction as shown schematically in Fig. [T] 
The corresponding vector potential in polar coordinates 
is Ar = 0,Ag = -JM. 

Let us assume, for simplicity, that without interaction 
with magnetization the electron spectrum corresponds 
to that of isotropic (in plane) Dirac cone. Then, the 
dynamics of the low energy states of the system in the 
presence of the vortex is described in polar coordinates 
by the Hamiltonian 
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where we have put h — 1. From now on and since we 
will not deal with interactions which can renormalize the 
Fermi velocity [53] we will put it just to unity: t; = 1. 
The usual ansatz 

vE'i(r,0) = e''ViM , vI/2(r,0)=e'('+i)V2(r), (2) 

allows to write the Schrodinger equation _ff 'f ~ E'i> with 
= (1'i,«'2)^ as 
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which can be reduced to 
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for both components. 

Equation ([4| is formally identical to the Schrodinger 
equation of a planar hydrogen atom with a nucleus of 
charge Ze^\^: 



A — 2m I (p 
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with the substitution e ^ E^ ~ {JM)^, 2mZe^ = 
JM(2l + 1). From this formal analogy we can write the 
answer for the eigenvalues of the energy: 



E'^irirJ) = {JMf 
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with rij. = 0, 1, 2... 

To ease the readout of the spectrum wc show a plot of 
the numerical solution of eq. ([!]) in Fig. [2] (center). The 
results reproduce the analytical spectrum given in eq. ^ 
but only for one sign of the angular momentum. There 
are no bound states for the other sign. For comparison we 
also plot the numerical results for the spectrum of a disc 
as a function of the angular momentum in the presence 
of a constant magnetic field (Fig. [2] (a)). 

We can see that for 2Z -I- 1 > (assuming JM > 0) 
there are localized states in the energy interval \E\ < JM 
with accumulation points at the boundary values E = 
zLJM. The numerical solution shows also that the wave 
functions for the localized states decay exponentially at 
large distances. For any value of I there is a zero-energy 
mode with rir — 0. 

For 2/ + 1 < we have the effective Schrodinger equa- 
tion for the Coulomb center with repulsion, and no local- 
ized states. 
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FIG. 2: Top: Energy levels as function of angular momentum 
at the surface of a topological insulator in a constant in-plane 
magnetic field. Center: Energy levels in the presence of a 
magnetic vortex. Bottom: Energy levels in the presence of a 
skyrmion configuration. 



The natural spatial scale in the problem is ^ 



\J\M 



which plays the same role as the correlation length for 
the Abrikosov vortices. 



Skyrmion configuration 

A novel situation arises when the (3D) magnetic coat- 
ing is in a skyrmion configuration |26j . 2D skyrmion 
patterns can exist when the magnetic anisotropy is weak, 
and they have been recently observed with Lorentz trans- 
mission electron microscopy in a thin film of Feo.5Coo.5Si 
[27j . The magnetization vector evolves from pointing 
along a given direction at the center of the defect to the 
opposite direction at large distances. The absolute mag- 
netization does not change, and at intermediate distances 
it has an in-plane component. A possible skyrmion tex- 



ture is 

M(r, e) = (sin /(r) cos 9, sin f{r) sin 6*, cos /(r)), (7) 

where /(r) can be chosen as a constant for simplicity. 

A numerical solution of the spectrum induced by a 
skyrmion texture is shown in Fig. [2] (bottom) . The cal- 
culation solves an effective one dimensional differential 
equation for each value of the angular momentum j28j . 
The figure also shows spectra for a constant magnetic 
field (top), and for the gauge field associated to a vor- 
tex (center). The skyrmion opens a gap in the spectrum 
and does not support a zero energy solution. Note that a 
zero mode appears when the in-plane component of the 
magnetization is subtracted. 



CONCLUSIONS AND DISCUSSION 

This work originated on the question of whether an 
effective gauge field would be induced by strain on the 
surface of a topological insulator. As it is known this is 
the case in graphene which can be seen in many respects 
of a precursor of the actual topological insulators. It 
is known that elastic deformations and geometrical cor- 
rugations in graphene can be described by and induced 
fictitious gauge field coupling to the electronic degrees of 
freedom The apparent paradox that the strains do 
not break time reversal symmetry while magnetic fields 
do is solved in the graphene case by the fact that the in- 
duced fictitious field couples to the two Dirac cones with 
opposite signs and time reversal symmetry is preserved in 
the complete system. This mechanism will probably oc- 
cur also in the case of weak topological insulators with an 
even number of Dirac cones at the surface. The response 
of a strong topological insulator to elastic deformations 
has to be different and remains to be studied. In this 
work we have seen an alternative way to generate gauge 
fields in a strong topological insulator and we have stud- 
ied the modification of the energy spectrum caused by 
the gauge fields. 

Magnetic coating of the surface three dimensional TI 
is often invoked to describe specific physical situations 
where a gap is needed in the 2D system [29H3T]. The 
discussion presented in this work shows that in specific 
situations the magnetic coating coating can give rise to 
a strong reorganization of the spectrum that has to be 
taken into account for each physical situation. 
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